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The effect of the charge fluctuations on the electronic and Raman spectrum of high 
temperature superconductors is examined, using the slave boson approach to the large 



O ■ Coulomb repulsion Ud on the copper site. Instead of the saddle point approximation 
(N 

^ . (6) 7^ for the slave boson, characteristic for various N ^ oo approaches, we confine 
Q^ - 

ourselves to the non-crossing approximation (NCA) diagrams for the Green func- 

^ : 

S ■ tions. In this way the effects of charge fluctuations and of the constraint of no double 
^ ' occupancy on the copper site on the shape of the electronic spectrum are studied 

o ■ 

^ ' primarily, while the slave boson diagrams responsible for copper-copper spin corre- 
^ i lations are intentionally not included. The novel feature of the charge fluctuation 
dynamics shows up in the slave boson spectrum as the plateau extending over the 
range of leV below u = 0. This is further reflected in the electronic Raman spectrum 
that we calculate for various combinations of the polarization of the incoming and 
scattered light. The Raman spectrum shows the characteristic featureless behaviour 
up to frequencies of the order of leV while the polarization dependence is also in 
qualitative agreement with experiment. 
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1 Introduction 



The attention which the copper-oxide superconductors have been attracting in recent 
years is due to their extraordinarily high transition temperature to the superconduct- 
ing state. The interest in normal state properties stems from the belief that a better 
understanding of these properties could give a clue to the superconducting mecha- 
nism itself. The Raman spectrum, and especially its continuous, almost featureless, 
background spreading up to the frequencies of the order of leV, is one of the most 
remarkable of these properties [jl], 0]. Considering the energy scales involved it is 
natural to ascribe this background to the electronic excitations in the system. In this 
work we propose a microscopic theory of this phenomenon, which is believed to be a 
common feature of all copper-oxide superconductors. 

Several possibilities for the origin of the Raman background were discussed pre- 
viously in the literature. We comment briefly on some characteristic ones. First, 
there are theories in which the Raman background is attributed to particular band 
crossing effects [Q. However, the one-dimensional bands used in the particular case 
of YBa2Cu307_j, (where they come from the CuO chains) are not expected to be a 
feature of the copper-oxide superconductors in general. Second, the "Nested Fermi 
Surface" (NFS) approach Q starts from the assumption that parts of the Fermi sur- 
face have a pronounced nesting property and proceeds by considering the effect of a 
weak Coulomb interaction. While this approach is certainly an interesting one, vari- 
ous Fermi surface experiments and band calculations indicate that the Fermi surface 
geometry varies significantly among the HTSC materials, which show the featureless 
Raman background. Third, the possibility was also considered that the magnetic ex- 
citations (paramagnons, related to the antiferromagnetic phase of the pure, undoped 
compound and common to most HTSC materials) are related to the Raman back- 
ground, for instance through the coupling of light to the two-magnon excitations. 
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However, the magnetic exchange constant J, measured by the polarized neutron 
scattering in pure materials, is not big enough (e.g. J ~ ITOO-R' in YBa2Cu306, 
and J ~ MOOfT in La2Cu04 0) to account for the excitations of the order of 
8000cm~^. Furthermore, the magnon (paramagnon) excitations presumably become 
more damped as the doping level is increased and the effective J is significantly 
suppressed 0, which has been considered and qualitatively understood by several 
theoretical approaches. On the other hand, it seems that the magnetic excitations 
indeed contribute to the Raman background around 4000cm~^. This contribution is 
pronounced [0, |§, ^ in pure and slightly doped materials and fades out on further 
doping. In fact it seems that this contribution is superimposed on an almost uniform, 
featureless background extending up to ~ leV, present in superconductors. 

Therefore, in order to understand the origin of this background, we turn towards 
the role of charge fluctuations and the effect of a strong Coulomb interaction, present 
in most models of HTSC's. The possibility that the charge fluctuations are responsible 
for the special properties of the normal phase was suggested by several authors from 
the very beginning. Varma and collaborators were seeking for the microscopic origin of 
their marginal Fermi liquid phenomenology in this domain. However, their approach 
was essentially a weak coupling one while most of the model parameter evaluations 
give an on-site repulsion Ud on the copper atom which is large relative to the copper- 
oxygen hybridization, Ud ^ ^o- Therefore, we use a strong coupling approach to 
calculate the effect of large Ud on the Raman spectra. However, the interatomic 
Coulomb repulsion term Vpd-, important in their theory, is not included in the present 
calculation. The effect of the local constraint imposed by the Coulomb repulsion on 
the Raman spectrum is essential in our approach. In that respect it shares some 



major physical points with the recent work by Bang [|Tl|] in which the t-J model is 
considered within the more usual approach of the saddle point approximation plus 
some fluctuation corrections in the slave bosons. 



2 Model and method 



Following many previous works on the copper-oxide superconductors (refs. |TJ, |T3|, |14 



and many others) we treat the large Coulomb repulsion Ua in the f/^ — oo limit using 



the slave boson formalism [15, 16, 17|. The application of this formalism is usually 
accompanied by the immediate use of the saddle point approximation (b) ^ for 
the slave boson field b. This approach has successfully accounted for some electronic 
properties of the system such as the in-gap "band" containing the Fermi level, the 
large electronic mass and the electronic dispersion of the form found in ARPES ex- 
periments. Furthermore, the integration of the high energy (few eV) fluctuations 
(cj ~ — 6d) of the boson field around the saddle point adds one more important 
ingredient to the physics of narrow band fermions - the copper-copper magnetic in- 
teraction J^d [111]- However, it is the intermediate frequency range, u < leV, that 
must be primarily considered for the explanation of the Raman background. In this 
region of frequencies, the saddle point approximation does not seem to be an advan- 
tageous starting point. In particular, as with any system of noninteracting fermions, 
the system of effective fermions that emerges on the saddle point level has no intra- 
band k = excitations, so that self-energy corrections for fermion lines and vertex 
corrections in the Raman scattering diagrams are required in order to get a finite 
Raman response for the frequencies of interest. Furthermore, the spectrum of the 
Gaussian fiuctuations of the slave boson field around their saddle point value shows 



[p!9| no sign of the continuum present in the Raman scattering. Also, some funda- 



mental difficulties of the saddle point approach |^0| are accompanied by some more 



practical ones. For example, the fc-dependence of the boson Green function and the 
corresponding /c-space integrals considerably complicate the calculation beyond the 
saddle point approximation requiring extensive numerical calculations. Actually, 
the intersite correlation of the bose field should be absent due to the local gauge in- 
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variance of the model in the slave boson representation, and to the slow fluctuations 
of the boson phase Therefore, we do not use the saddle point approximation for 
the boson field b in our calculations. Instead, we consider the high order diagrams in 
the boson field, in which the electronic density of states, similar in shape to the one 
obtained through the saddle point approximation (with the aforementioned in-gap 



resonance), and the electronic damping, emerge on equal footing E^. This approach 



is a direct extension of the treatment of Kroha et al. of the Anderson impurity 
problem. The model that we use for holes in Cu02 planes is the well known p-d 
model P3, 



H = Ho + Hhyb + Hd , (1) 

Ho = ^Up{clfi^CxHa + c\na'^yna) + (^dd^^d^a} , (2) 
n,(T 

Hhyb ^ ^, to\dfj„(^Cxna ~l~ Cyfl^ Cxn—l^a (^yft—lyo) ~l~ h.C.^ , (3) 

n,(T 

Hd = Yl Uddl^dft^di^d^i , (4) 

in which, apart from the Coulomb repulsion Ud, the parameters e^, ep and to denote 
respectively the energies of the hole in Cu Sd^^-yi and O 2p(cr) orbitals and the copper- 
oxygen hybridization. Transformed to the fc-space, Hhyb, acquires the following form, 
with the lattice constant set to unity: 

Hhyb = T.{dlhuJ.k + CykJykl + iikXk + ik/ykl'^k'r} ' (5) 

k,cr 

t.k^to{l-e'^'^), t^^, ^ toil - e-^'y) . (6) 

The Hamiltonian simplifies on introducing the bonding and non-bonding combina- 
tions of the hole operators on oxygen sites, 

_ _ ^xk^xk + ^yk^yk _ ^yk^xk ~ ^xlfyk 

Pk-^bk- /I . 12 , I =^ ' ^nfc - fiYZ\2 , I + . 12 ' 
y \ ''xk \ ^ \ ''yk \ y \ ''xk \ ^ \ ''yk \ 

Only the bonding combination has an overlap with the copper orbitals. In terms 
of these operators, the part of the Hamiltonian which describes noninteracting holes 



becomes 



Therefore, from now on we omit the non-bonding states from our discussion [p5| . 



In the slave boson representation |T5[ in the Ud ^ oo hmit, the fermion operator 
fl^ is introduced to describe a singly occupied, and the boson operator bt is introduced 
to describe an empty orbital state at the copper site n. The substitution dfia — >■ 
bfifna, which leads to the Hamiltonian in the new representation, is accompanied by 
a constraint, which separates the physical states from many non-physical ones in the 
enlarged Hilbert space. In terms of the number operator Qn, 

Qn=E/nl/n- + fe, (10) 

a 

which commutes with the Hamiltonian, this constraint is written as 

Qn = l. (11) 

In the path integral formulation, the constraint is introduced through the "Lagrange 



multiplier" A^j in the partition function El 



Z{Qn = 1) = <j n y ^^/^ [ Z({A})e^SH ^ (12) 



Z{{\}) = J Vb^VbVf^VfVp^Vp e-io''^W'^" , (13) 

d d d 

^(r) = E/i(^ + ed)/na + EpL(^ + ep)P- + E4^&n (14) 

n,cr n,(j n 

+H,y, + E A„^(E flfri. + bibn) . 

The saddle point approximation for the integration over in (|T2p fixes A^ to the 
value A which is, due to the translational symmetry of the system, equal on all copper 
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sites. Within this approximation, the constraint ([TT|) is satisfied on the average, 

m = -\^^ Hzm)) i(A,=A)= 1 . (15) 



Following the approach which was previously tested by Kroha et al. I^S] on the 
[/ — s> oo Anderson model, with the result that it is in good agreement with those 
obtained by other methods, we neglect the fiuctuations of around A. This means 
that the states with different from unity become allowed. However, the effective 
Hamiltonian, 

' t(k) ^''^ 
+7= E[pL E blfu,. + E fUMk.]} , 

with Xfi fixed to A, conserves Qn, so the states with different Qn are not dynamically 
mixed (A^ stays for the number of unit cells in the crystal). Therefore, we expect 
that fixing (Qn) = 1 in our calculation (which itself conserves the Q-charge) will 
predominantly extract the dynamics of the system with Qn = 1- 

It can be easily checked that, in the absence of the 6-condensate, the one-particle 
Green functions for b and / fields are site-diagonal (or constant in the fc-space) to 

all orders in the perturbation {1/VN) 'Zk,^,Ji^)ip\b\h+qa + fl+^J^qPka) in (HI)- 
simple words, it is impossible to create (or annihilate) a 6 or / particle at some site, 
and then propagate it to some other site, since that would violate the conservation 
of Q^. 

Further following Kroha et al. we confine ourselves to a non-crossing approx- 
imation (NCA) for which the free energy is represented by diagrams given in Fig. |I[ 
Fig. ^ shows the corresponding diagrams for the full Green functions for 6, / and p 
fields. The set of Dyson equations reads as 

g,{uj)-^ = Gl{uj)-^ - S,{uj) , (17) 
->o 



gp{u,k)-' = gl{uj)-^-s,{uo,k) , (19) 

with self energies E^)^ Sp and Sj of the following form 

1 - ^ - 2 

Sbi,^) = ^^\t{k) ab{uj,k) , ab{uj,k) = -Y,Qp{^,k)gf{uj + u) , (20) 



I!p{u,k) =\t{k) \^ ap{uj) , ap(u;) = ^ ^ ^^(1^)^/(0; + z/) , (21) 



^/(^) = ^ E I m |2 a/(^, fc) , af{uj, ^) = ^ E - ^) • (22) 

The bare Green functions are 

Q'A^) = -. — — — , ^.°M = -^, ^JH = . — 

^cj — A — + /i zcj — A ^ zcj — ep + /i 

where /i denotes the chemical potential for holes. The system of equations (|l^-(^) 
has to be solved self-consistently keeping at the same time {Qn) = 1- The number of 
holes per unit cell is fixed by the doping level of the system, Up + Ud = I + S. 

In solving the system of equations (p!7| ) -(p2| ) , we transfer them, by means of usual 
methods of analytic continuation |^ to the expressions in which the integrals over 
various spectral functions appear. In the numerical integration over frequency a mesh 
of thousand points per eV is used. Starting with bare Green functions we iterate the 
complete set of equations until the result does not change any more. Depending on 
the iteration scheme, the rate of convergence may change considerably, but the same 
results emerge at the end. For numerical purposes, we use narrow Lorentzians instead 
of the delta functions for the spectrum of the bare Green functions for fermions. Their 
widths are typically several times of the mesh size. Once the Green function Qp{uj, k) 
is known we obtain the Green function for the real hole on the copper site by using 
the relation 

g^iu, k)-' = cjpiuj)-'- I t{k) |2 gl{u) , (23) 

shown diagrammatically in Fig.^. 
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3 Hole spectrum and the local density of states 

In this section we present the results for the Green functions, which are obtained on 
solving the system of equations (]T7|)-(p3D. For further discussion, we choose parame- 
ters Apd = €p — ed = 4.0to (^0 is of the order of leV), 5 = 0.13 extra holes per unit 
cell, and set ep = for convenience. The values for A and fi that correspond |28| to 



this choice of parameters are A = 3.75to and /i = — l.lSSto- The chemical potential is 
positioned between the energies ed and of the copper and oxygen orbitals and the 
temperature is set to 0.0135to- 

For the /-fermion spectrum we get a simple Lorentzian-like peak (Fig. P, with 
about 80% of spectral weight in it, positioned on the Fermi level and with the oc- 
cupancy nf = 0.8 (both spin projections). The structure of the 6-boson spectrum is 
much richer, with a high frequency peak centered at ~ A, and with a tail extending 
to higher frequencies. For lower frequencies, we find a.) the asymmetric peak at 
w ~ and b.) the broad plateau starting from u ^ and extending beyond -leV. 
These features give rise to the following physical hole spectrum. 

The hole density of states per spin (hereafter DOS), — (I/A^tt) IniG'^(A;, tu), 
accumulates (see Fig. |) around energies ( "p-band"), ed ( "d-hand") and /i ( "in-gap- 
hand"). While the accumulation near ep and occurs for physically obvious reasons, 
the third region may be related to the band that occurs around the renormalized 
copper energy level in the saddle point approximation |T2[ for the h field. Experimen- 
tally, the "in-gap" states show directly as a bump in the PES/IPES spectra upon 



doping. The dispersion that occurs in the "band" around is characterized by the 
shift of the peaks in \m.Gp{k, u) towards higher hole energies along the TX line in the 
Brillouin zone. This is in qualitative accordance with a simple tight-binding calcula- 
tion for the f/rf = case, except for the fact that now the peaks are of considerable 
width. The scattering of fermions on the slow boson field resembles the scattering on 
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the random potential and, similarly to the crystal disorder, increases the relaxation 
rate for the fermions. The dispersions of the other two "bands" go in the opposite 
direction (see Fig. |^). Relative intensities of DOS for the oxygen and copper site 
show that "bands" close to and fi are dominantly of oxygen character, while the 
"band" around is, as expected, mainly of copper character. The overall result is 
qualitatively similar to the DOS obtained by carrying out the Gaussian corrections 
beyond the saddle point approximation in the 1/N expansion for the generalized 
A^-component model |^ . 



At this point the total number of states on the copper and the oxygen atom 
can be found. The total number of states per spin at the oxygen "site", given by 
the integral of — (l/7r)ImG^(fc, cj), is 1, as it should be. However, the number of 
states per spin per copper "site" should be 0.5, as one state per copper site is shifted 
to infinity in the Ud —>■ oo limit, while the integral over ImG^ gives a somewhat 
larger number, 0.5 • [1 + (1 — Uf)]. This deviation can be readily traced back to the 



= const approximation which we use and to the decoupling |^ of b and / fields 
that takes place in the calculation of dp. Indeed, after the analytic continuation to 
real frequencies, the equation (|2T| ) gives 

-/^Imap^M = - J ^{lmG^{u)cth{^) -lmGfiu)th{^)} 

(24) 

1 I 1 I i-fif 

= 2^/ + ^fe = 2 + 2 ' 

These extra (1 — nj)/2 states enter into the DOS through the hybridization equations 
(|T9|) and (p3D producing the aforementioned deviation. 

One feature of our results which seems alarming when compared with the ARPES 
measurement of the Fermi surface geometry is that the chemical potential /i, as ob- 
tained by the presented selfconsistent calculation, is positioned above the peaks in 
lmG{k,uj). It means that there is practically no Fermi surface in our calculation, 
although /i lays inside the in-gap band since the tails in lmG{k,uj) extend above /x. 
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The reason for this may be traced back to the position of the peak in lmGf,{uj) which 
is close to = 0. The precise position of the peak determines the Fermi surface 
geometry and the shift of this peak of the order of O.OSeV to lower frequency would 
cure the problem. We believe that the calculation beyond NCA may produce such a 
shift restoring the Fermi surface consistent with the experiment. On the other hand, 
at the energy scales larger than a tenth of an eV the major features of the spec- 
tra would presumably remain unchanged in the improved calculation. These energy 
scales are of the prime interest in further discussion. This concerns in particular the 
plateau in the boson excitation spectrum which extends over the range of leV. 



In the calculation of the Raman spectra this plateau implies the occurrence of the 
flat background. 

4 Raman spectrum 

Before describing the calculation of the Raman spectrum we wish to recall the simple 
results obtained with the electronic relaxation rate introduced artificially |3^. Using 
the fermions with the Fermi liquid like relaxation rate 'y ^ u'^/T gives a Raman spec- 
trum with a maximum at the frequency of the order F. Further on, such spectra with 
a maximum are obtained for a wider class of the non-Fermi liquid electronic relax- 
ation rate laws 7 ~ r{uj/r)°' (with a > 0, a 7^ 1; F always being the characteristic 
energy scale involved) This is in contradiction with the experimentally observed Ra- 
man scattering intensity which does not diminish (or which even increases) towards 
high transfer frequencies of leV. The phenomenological marginal Fermi liquid ansatz 
of Varma et al. [0] (a = 1 in the above 7) has no such parameter F and a linear 
damping rate for electrons leads to the constant Raman background. 

In what follows we do not use any artificial assumption for the electronic spectrum. 
Instead, we use the boson and the fermion propagators obtained by the numerical 
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calculation of the preceding section. While it is difficult to draw the correspondence 
between the aforementioned phenomenological calculation and the ours, it is inter- 
esting to note that our plateau in the slave boson spectral function is the feature 
which leads to the nearly constant Raman background similar to that obtained in the 



marginal Fermi liquid phenomenology 



The calculation of the Raman spectrum proceeds as follows. The coupling of the 
electromagnetic field to electrons is introduced, as usual, by the substitution 

in the interatomic hybridization term of the original Hamiltonian. This is followed by 
the expansion of the exponentials up to the terms quadratic in the vector potential A. 
The terms p ■ A and nA"^ describe the coupling of the electromagnetic field to the hole 
momentum and to the bond charge. As usual, these terms give rise to three different 
contributions Tq, Th and Tc to the transition matrix for the Raman scattering. They 
are shown in Fig. |^ and correspond to 

a) nA"^ scattering; 

b) {p ■ AY scattering where electron first absorbs the incoming photon with fre- 
quency Ui and then emits the outcoming photon with frequency 1^5; 

c) similar to b) but with absorption and emission processes taking place in the 
reverse time order. 

The probability for the transition between the initial state |0) (for simplicity we 
use zero temperature formulae to illustrate further steps) and the final states |/) 
consistent with the Raman processes described above is equal to 

= y E K^a + + Tc)of I' S{Eo - Ef) . (26) 
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This gives six (three diagonal and three nondiagonal) contributions to the transition 
probabihty, 

W = Waa + W,, + W,, + Wab + Wac + W,, , (27) 

W,^ = {2-5i,)lY.(T:Tj+T;T,)5{E,-Ef) , z,j = a,b,c. (28) 

The Raman spectra that emerge after the numerical integration are shown in 
Fig. p. The featureless behaviour that characterizes the resulting Raman spectrum 
extends approximately up to the frequency of the separation of the "in-gap-band" and 
the "p-band". This separation is of the order of leV. This agrees with experiment, 
while scatterings for larger frequencies are not experimentally probed. Above that 
frequency, our calculations give a sudden increase of the Raman intensity, due to the 
"interband" processes. Coming back to the spectrum below leV, we note that the 
intensity of the Raman spectrum changes significantly as a function of the polarization 
of the incident and the scattered light. This is due to the polarization dependence 
of the electron- light interaction which comes from the nearest-neighbour hopping 
approximation (see the Appendix for some details). An example of the relative size 
of contributions which come from various terms in (^) is shown in Fig. |^. The relative 
importance of various terms changes with polarization (however, the diagonal terms 
almost always dominate over the nondiagonal ones). In particular, the scattering via 
the nA^ coupling does not contribute at all for the polarization combination {x',y'). 
In general, the polarization dependence of the intensity of the Raman spectrum shows 
good qualitative agreement with experimental results. The minor difference in that 
respect is that the intensity of our spectra for (x, x) and (?/', y') polarizations are 
somewhat lower that the experimental ones This may point to the necessity 
of including direct oxygen-oxygen hybridization terms (as well as the corresponding 
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electron- light coupling terms) in the initial p-d model. 



In summary, we have calculated the hole spectrum in the Cu02 plane, using the 
slave boson formalism to implement the constraint of no double occupancy on the 
copper site, and without resorting to the usual saddle point approximation for the 
slave boson field. A three band density of states structure is obtained, with the 
Fermi level lying in the "in-gap" band. The corresponding energy scales show up in 
the Raman spectra, featureless in frequency up to the value of the gap between the 
"in-gap-band" and the "p-band" . The measured polarization dependence agrees well 
with our calculation, which is based on the assumption that the Cu-0 hopping to is 
the dominant one in the Cu02 planes. 



This work was supported by the Croatian Government under project 1-03-275 and 
by the EU-Croatia project CI1*0568-C(EDB). 



A Appendix 

The polarization dependence (and the corresponding symmetry properties) of the 
Raman spectrum result directly from the polarization dependence of the electron- 

— * 

light coupling terms in the fc-space representation. The coupling constants g\ and g2 
correspond respectively to p- A and nA^ terms. In the p-d model in which only Cu-0 
hopping is included these coupling constants have the following form: 

^^'^ / I f . 12 , I f . 12 ^ ' ^^^> 

Y I ''life I "I" I Sfe I 

<k ^ /I.. 12 , I.. 12 — ' (30) 



yk I 

dn [(ej • x){es ■ x) - jei ■ y){es ■ y)]t:ckKk 
07. t oc . — , (31) 

xk ' ' yk \ 
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dp [i^i ■ ^)i^s ■ x) I P +(e/ ■ y){ei ■ y) \ t ^ p] 
k ^ , , (32) 

V I ^xk P + I ^yk P 

where 

^.k ^ + e-''^) , V = + e-''^) . (33) 

Here e} and es denote polarization vectors of the incident and scattered hght (in 
gi expressions ep=e} or e^). The superscripts dp and dn stand for the case of the 
couphng of the hght to the copper-bonding and the copper-nonbonding oxygen terms 
respectively. 
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Fig. 1: First few diagrams for the free energy within the NCA. 

Fig. 2: The diagrams for the Dyson equations for 6, / and p fields within the NCA 
approximation. / and h are site-diagonaL Bold lines denote the full Green 
functions. 

Fig. 3: The diagram for the Green function for the hole on the copper site. 

Fig. 4: The spectra for the (a) /-fermion and (b) 6-boson. 

Fig. 5: DOS for (a) the oxygen site and (b) the copper site holes. 

Fig. 6: In-gap band dispersion (the one containing the Fermi level) along the VX line. 
The spectra for (a) the oxygen site and (b) the copper site are shown. 

Fig. 7: Diagrams for the lowest order contributions to the T-matrix for the Raman 
scattering. 

Fig. 8: Transition probability for the Raman scattering for various polarizations of the 
incident and scattered light as the function of the transfer energy uj. 

Fig. 9: The contributions {Waa, Wbb, Wcc in (a) and Wab, Wac Wbc in (b)), coming from 
various T-matrix terms, to the total transition probability W are shown for the 
{x, x) polarization. 
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